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We consider the issue of chiral ordering in spin chains for an arbitrary value of the spins 
S. By use of bosonization according to a scheme developped by H.-J. Schulz, we obtain the 
phase diagram of the chain with XY exchange couplings between nearest and next-to-nearest 
neighbors. We obtain a satisfactory picture including a so-called chiral spin nematic phase 
which is gapless, has long-range chiral order and incommensurate spin correlations. We 
perform a stability analysis of this phase and point out that this analysis is in conflict with 
existing DMRG results that shows a difference between integer and half-integer spin case. 

§1. Introduction 

Quantum antiferromagnetic spin chains display a variety of phases that have 
no classical counterpart. This variety is even increased if we consider the effect of 
frustration. Many studies!]) have been devoted in the past to the simple AF chain 
with nearest neighbor exchange J\ and next-to-nearest neighbor J2. In the isotropic 
case, the Hamiltonian is then simply : 

H = J!J2 (SnS n+ l) + J2J2 ( S nSn+2) , (M) 

n n 

where = ± iSf[ is a spin operator at site n and there are competing antifer- 
romagnetic interactions Ji, J2 > which introduces frustration in the model. In the 
spin-1/2 case, for small J2, there is a spin-fluid phase whose effective theory is that 
of a massless free boson. It has quasi-long range spin order with algebraic decay of 
spn correlations. For larger values of J2, the ground state is spontaneously dimer- 
ized : this appears through a quantum phase transition of Kosterlitz-Thouless type 
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and eventually incommensurability develops within this gapped phase for even larger 
values of J2 but without any additional phase transition. The situation is remarkably 
different if we now consider XY exchange : 

% = Jl 5^ (^n^n+l + ^n^n+l) + <^2 + ' O-'ty 

n n 

Starting from the limit with large J2, Nersesyan et al.9 used a mean- field treatment 
and then bosonization to predict the occurence of a new phase with many uncon- 
ventional characteristics. In this limit, they predicted that there is long-range chiral 
order : 

((S n A S n +i) z ) / 0. (1-3) 

There are also local spin currents polarized along the anisotropy z-axis. This phase 
is gapless and there are incommensurate spin correlations that decay algebraically 
with an exponent which they found to be 1/4. The existence of this phase has been 
recently demonstrated numerically. Such a phase has been disclosed in a ladder 
model formulated first as an array of Josephson junctions which is equivalent to a 
spin- 1/2 half model &. this ladder has square plaquettes so it is not frustrated as 
model Eq. (|l-2[) but frustration is introduced by half a flux quantum piercing the 
plaquettes. Then a study & using the DMRG algorithm gave evidence for this phase 
in model Eq.(|l-2j) : the spin fluid phase is stable up to J1/J2 ~ 0.33 then the chain 
undergoes dimerization and at J1/J2 ~ 1-26 there is a second transition to the chiral 
critical phase. This new phase has also been reportedH^ in the S=l chain with the 
same Hamiltonian Eq.( |l-2j ). Here the situation is even more richer. When J2 = 0, 
we are in the XY phase which destroyed immediately by adding even an infinitesimal 
J2, the phase that appears then is the celebrated gapped Haldane phase. This phase 
resists the perturbing influence of J2 for a while but at J1/J2 ~ 0.47 there is a 
phase transition to chiral order but the gap remains nonzero. Then very close, at 
J1/J2 ~ 0.49, there is a distinct transition to a critical phase with chiral order as in 
the S=l/2 case. So the integer S=l case has an additional phase w.r.t. the S=l/2 
case, a gapped chiral phase. This difference persists for higher spinsB. For S=3/2, 
there is an intermediate dimerized phase which is replaced in a single tranistion by 
by chiral critical phase. For S=2, the Haldane phase is destroyed by two successive 
transitions as for S=l. To understand this pattern of phase transition we use& 
the bosonization technique which has been adapted to the case of generic spin-S 
by SchulzB*. This method is able to capture the phase diagrams as a function of 
exchange anisotropy as well as single-ion anisotropy and it correctly captures the 
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difference between integer and half-integer spins. 

§2. Weak coupling limit 

The idea is to write down each spin-S as a sum spins 1/2 : 

25' 

Sn=X>ir (2-1) 
<x=l 

Each spin 1/2 is then bosonized by the standard technique. There is a bosonic field 
4> a , a = 1 . . . 2S associated with each spin. The first possibility is to treat J<i as a 
perturbation. So we start from an isolated spin-S chain and write its effective theory 
in terms of the bosons <f> a . Due to the appearance of couplings sZ, there are 
operators that induce gaps for some linear combination of the basic bose field. More 
precisely, only the "acoustic" mode remains massless : 

* = -i (fi + - + <P2S) ■ (2-2) 



The effective Hamiltonian for the acoustic mode is thus a simple free theory : 

Hxyi ~\( n2 + {dx$?) , (2-3) 

where 77 is the canonical momentum conjugate to The coefficient v is an unim- 
portant velocity and we have used the conventional name "XY1" coming from the 
standard S=l chain phase diagram. The spin operator can be expressed in terms of 
the field which is dual to ^ : 



S ± ~ (-lf /a exp {±i^fir/2S <9j . (2-4) 

This expression shows easily that the XY spin correlations decay algebraically with 
an exponent r\ = 1/AS. This is in agrrement with numerical findings 0>. In the S=l 
case it is also in agreement with work by Kitazawa et al. If we now bosonize 
the perturbation J2, we find that there is a simple renormalization of the previous 
free hamiltonian Eq.( |2-3| ) but in addition vertex operators appear in perturbation 
expansion. For integer spin, the most relevant operator appear at S th order and it 
appears at 2S th order for half-integer spins. The effective theory is then : 

H~ V - (KIT 2 + 1 (d^) 2 ) - ^ cos (P$) , (2-5) 



where f3 = V8ttS for integer spins and (3 = V32irS for half-integer spins. The 
Luttinger parameter K is obtained in perturbation K = 1 — (A/tt)^/ J\ + 0(Jf )• 
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This effective leads then immediately to the phase diagram of the spin-S XY chain 
for small J2/ J\- The scaling dimension of the vertex operator in the effective theory 
Eq. (|2-5| ) is K(3 2 /4tt and thus is irrelevant for small J2, the XY1 phase will have thus 
a finite extent. With increasing J2, the vertex operator becomes relevant and drives 
the system towards a massive phase through a KT transition. Depending upon the 
spin parity, it will be the Haldane or the dimerized phase. If we take seriously 
the approximate formula for K, we deduce J2/J1 ~ 0.29 at the KT transition for 
S=l/2, which compares quite favorably to the numerical estimate of ~ 0.324. We 
also predict that S=l is special : in this case indeed the operator is marginal for 
J2 — * hence the instability takes place immediately upon switching an infinitesimal 
value of J2, as seen in DMRG studies. 



§3. Zigzag limit 



We now turn to the opposite limit J2 3> J\. Then we have a two-leg spin-S XY 
ladder coupled in a zigzag way. We first bosonize the two independent legs when 
J\ = 0. So each of the chains can be treated as in the previous section. There are now 
two acoustic modes $1 and <&2 that are the effective low-energy degrees of freedom. 
It is convenient to introduce the symmetric and antisymmetric combinations of these 
two modes : 

$± = ■^=(^±#3). (3-1) 
The leading contribution to spin correlations comes from the following operator : 

St * (-ir /a exp (±i^f/2S a ) , (3-2) 

where a = 1, 2 and Q a are the fields dual to <P a . When J\ = the two fields <P a are 
free and massless. Introducing Ji, we find the effective theory : 

U - \ E { n a + (W 2 ) + 9 d x 9 + sin (^1 ©-) , (3-3) 

where Q± are fields dual to <I>± and II ± are canonical conjugate to <I>±, and g = 0{J\). 
The operator perturbing the free part in Eq.( |3-3| ) is a parity symmetry breaking 
term with nonzero conformal spin (=1). Its effect is thus highly nontrivial. A simple 
perturbation with nonzero conformal spin is given by the uniform component of the 
spin density d x <I>. In this case we know its effect : it induces incommensurability. 
We treat theory Eq.(p-3[) by following exactly the method of nersesyan et al. We 



Chiral Order in Spin Chains 



5 



just decouple : 

d x O+ sin (J^®-) ~^ K < d x@+ >+V< sin (^f -) > ■ ( 3 ' 4 ) 

We then impose self-consistency. The "+" sector remains massless while the "-" 
sector is massive due to the vertex operator in Eq. (]3-4j ). The results that follow 
are then close to the original findings for for S=l/2 critical chiral phase. The most 
notable difference is that spin correlations decay algebraically with a spin-dependent 
exponent : 

e iqx 

(Sl(x)S-(0))~-rr [ j m , a = 1,2, (3-5) 

with q — ir/a ~ (Ji/^) 45 ^ 45-1 )- The exponent 1/4 of the S=l/2 case is thus a 
special case of rj = 1/8S. This formula is in good agreement & with measurements 
by DMRG up to S=2. There are nontrivial spin currents in the ground state : 

(JL) = (JL) = ~vJ^ (d x O+) + 0. (3-6) 
In the language of spins, this means intrachain currents : 

((S a , n A S a , n+1 ) J ^ ^ (d x O + ) ^ 0, a = 1,2, (3-7) 
as well as interchain currents : 

Ji((Si,n /\ S 2 , n ) J oc (sin (J^&- W 0, (3-8) 
this also means long-range chiral order. 

§4. Stability analysis 

We certainly expect that when increasing the coupling J\ the "+" sector of 
theory fl3-3| ) will not remain massless forever because we know from the analysis of 
the previous approach in the small J2 limit that there is a massive phase. It is likely 
that some vertex operator will become relevant at some finite value of J\j 1 J%. We 
try to find this operator by analyzing the symmetry properties of the theory. Taking 
into account the symmetries of the mean-field Hamiltonian, we find the following 
effective theory : 

n+ ~ ~ {xnl + 1 (a^ + ) 2 ) + Kd x e+ - ^ cos ( 7 <p+) , (4-i) 



where 7 = V16ttS for integer S and 7 = V 64irS for half-integer spins. When this 
vertex operator becomes relevant, the system is gapped and still incommensurate. 
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This is consistent with the observation of a "chiral Haldane" phase in numerical 
studies for integer spins. For half-integer spin, this predicts a phase which is dimer- 
ized (due to the value at which <P + is pinned) and has chiral LRO. We are then 
forced to speculate that there is then an additional Ising transition at which chiral 
LRO disappear to make contact with the small J 2 limit. 

§5. Conclusions 

The bosonization approach is able to reproduce the phase diagram of the XY 
J\ — J2 chain for integer spins. For all spins it correctly predicts the existence of the 
chiral critical phase. The spin correlations decay with exponent 1/8S. For the half- 
integer case, this approach predicts a dimerized phase with incommensurability and 
chiral order which is apparently not seen. It is possible that this signals a failure of 
the mean-field decoupling. For integer spins, there is good agreement with a large-S 
study©. A possible candidate© for the chiral critical phase for S=l is CaV204 
which has J\ ~ J 2- This would require a large single-ion anisotropy© to escape 
from the double Haldane phase. 

Acknowledgements 

We would like to thank T. Hikihara and M. Kikuchi for information on their 
recent results. 

References 

[1] F. D. M. Haldane, Phys. Rev. B25, 4925 (1982); Erratum ibid. 26, 5257 (1982).; T. 

Tonegawa and I. Harada, J. Phys. Soc. Jpn. 56, 2153 (1987). 
[2] A. A. Nersesyan, A. O. Gogolin, and F. H. L. Essler, Phys. Rev. Lett. 81, 910 (1998). 
[3] Y. Nishiyama, Eur. Phys. J. B17, 295 (2000). 

[4] T. Hikihara, M. Kaburagi, and H. Kawamura, Can. J. Phys. 79, 1587 (2001). 
[5] M. Kaburagi, H. Kawamura, and T. Hikihara, J. Phys. Soc. Jpn. 68, 3185 (1999). 
[6] T. Hikihara, M. Kaburagi, H. Kawamura, and T. Tonegawa, J. Phys. Soc. Jpn. 69, 259 
(2000). 

[7] T. Hikihara, M. Kaburagi, and H. Kawamura, Phys. Rev. B63, 174430 (2001). 

[8] P. Lecheminant, Th. Jolicoeur and P. Azaria, Phys. Rev. B63, 174426 (2001). 

[9] H. J. Schulz, Phys. Rev. B34, 6372 (1986). 
[10] F. C. Alcaraz and A. Moreo, Phys. Rev. B46, 2896 (1992). 
[11] A. Kitazawa, K. Nomura, and K. Okamoto, Phys. Rev. Lett. 76, 4038 (1996). 
[12] A. K. Kolezhuk, Phys. Rev. B62, R6057 (2000). 
[13] H. Kikuchi, M. Chiba and T. Kubo, Can. J. Phys. 79, 1551 (2001). 
[14] T. Hikihara, Can. J. Phys. 79, 1593 (2001); J. Phys. Soc. Jpn. 71, 319 (2002). 



